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Abstract 



Models of disorder with a direction (constant imaginary vector-potential) 
are considered. These non-Hermitian models can appear as a result of com- 
putation for models of statistical physics using transfer matrix technique or 
describe non-equilibrium processes. Eigenenergies of non-Hermitian Hamil- 
tonians are not necessarily real and a joint probability density function of 
complex eigenvalues can characterize basic properties of the systems. This 
function is studied using the supersymmetry technique and a supermatrix o~- 
model is derived. The cr-model differs from already known by a new term. 
The zero-dimensional version of the cr-model turns out to be the same as that 
obtained recently for ensembles of random weakly non-Hermitian or asym- 
metric real matrices. Using a new parametrization for the supermatrix Q the 
density of complex eigenvalues is calculated in 0D for both the unitary and 
orthogonal ensembles. The function is drastically different in these two cases. 
It is everywhere smooth for the unitary ensemble but has a 5-functional con- 
tribution for the orthogonal one. This anomalous part means that a finite 
portion of eigenvalues remains real at any degree of the non-Hermiticity. All 
details of the calculations are presented. 
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I. INTRODUCTION 



Physics of disordered metals and semiconductors has been attracting a considerable 
attention during several decades. Various interesting phenomena were discovered experi- 
mentally and found a theoretical explanation. Rather simple models of a particle moving in 
a random potential can be used to describe such different effects as Anderson localization!, 
mesoscopic fluctuations!!, Integer Quantum Hall Effect!, and many others. 

Although the phenomena can occur already at a weak disorder, a simple perturbation 
theory in the disorder potential is not sufficient for their quantitative description. A proper 
theory is often based on summing certain classes of diagrams (cooperons and diffusons)!'! 
but in more complicated cases one has to use essentially non-perturbative methods like the 
supersymmetry technique! based on mapping of the disorder models onto a supermatrix a- 
model (for a recent review see Ref .! and references therein) . A disordered physical system can 
include a magnetic field, magnetic and spin-orbit impurities, etc. However, these additional 
interactions are included into the calculational schemes without considerable difficulties. 

By now, the diagrammatic expansions and the supersymmetry technique give a possi- 
bility of getting explicit results for most of the disorder problems. In addition, the super- 
symmetry method was applied for calculations with random matrices!, which resulted in 
application of the method in nuclear physics and quantum chaos where the random matrix 
theory (RMT) had been the basic computational tool (for a review see, e.g. Refs.ll2h0). 
Recently, a supermatrix cr-model was derived for ballistic billiards averaging over either rare 
impurities^! or energyB. So, the way of studying all these interesting problems appears 
quite clear, although in some cases one can encounter certain technical difficulties. 

The systems mentioned above are described by quantum mechanical Hermitian Hamil- 
tonians. After averaging over disorder the systems involved are invariant with respect to 
inversion of coordinates. Sometimes, in order to describe the decay width of eigenstates, 
non-Hermitian Hamiltonians are used. This approach is popular in study of quantum dots 
coupled to leads. Of course, the Hamiltonian of the whole system of the dot with the leads is 
Hermitian but it is often convenient to exclude the leads from the consideration by integrat- 
ing out degrees of freedom related to the leads. As a result of such an integration one comes 
to an effective non-Hermitian Hamiltonian of the dot containing imaginary energiesS. This 
type of the non-Hermiticity can be easily included into the scheme of the supersymmetry 
technique as well as into diagrammatic expansions and many results have been obtained 
explicitly!. 

In a recent publicationEl Hatano and Nelson considered another type of non-Hermitian 
Hamiltonians with a disorder, namely, Hamiltonians with a constant "imaginary vector 
potential". In other words, the Hamiltonians contain not only the second order derivative 
over space coordinate but also the first order derivative with a real coefficient. The model 
appears as a result of mapping of flux lines in a (d + 1) -dimensional superconductor to the 
world lines of d- dimensional bosons. Columnar defects produced experimentally by energetic 
heavy ion radiationll! in order to pin the flux lines lead to the random potential in the boson 
system, whereas the component of the magnetic field perpendicular to the defects results in 
the constant imaginary vector potential^. 

Already qualitative arguments!!! indicate that the presence of the imaginary vector po- 
tential can lead to new effects. In particular, a one-dimensional chain of the bosons has to 
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undergo a localization-delocalization transition; this result was also checked by a numerical 
computation. In "conventional" (without the first order derivative) disordered systems tran- 
sitions in one dimension do not occur and therefore the model with a direction belongs to a 
really new class of systems that have not been studied yet. It is argued that the localized 
states should have real eigenenergies whereas eigenenergies of the extended eigenstates may 
have a non-zero imaginary part. 

The importance of investigation of such systems becomes even more evident if one recalls 
that e.g. the equation for heat transfer with a convection has a term with the first order 
derivative. One can imagine a situation when quantum hopping of a particle from site to 
site of a lattice has a different probability depending on direction. The presence of the 
first order derivative in the Hamiltonian just corresponds to the introduction of a certain 
direction. The non-equivalence of the directions can be provided by coupling to another 
subsystem with broken inversion symmetry playing the role of a reservoir; this reservoir 
may be out of equilibrium. The classical analog of the disordered models with a direction 
(so called, directed percolation) has been discussed in the literature^. 

Another problem where one comes to a stochastic equation containing first order deriva- 
tives is the problem of turbulence in flow dynamics. It is generally believed that the most im- 
portant features of the turbulence can be described by the so called noisy Burgers equationil, 
which is a non-linear equation with a white noise random force. Besides its application in 
the flow dynamics this equation is used as a toy model by field theorists due to a striking 
analogy between the constant flux states in turbulence and some anomalies in quantum field 
theories^!. The Burgers equation is eauivalent to the Kardar-Parisi- Zhang equation intro- 
duced to describe the crystal growth". The non-linear Burgers equation can be reduced 
through a Hopf-Cole transformation to a linear (d + 1) dimensional equation with a random 
potential and time playing the role of the additional dimension. This equation has a first 
order time derivative and there have already been an attempt to solve it using the replica 
method@. The noisy Burgers equation can also be reduced to a quantum spin model with 
a non-Hermitian Hamiltonianl!. Recently, some interesting results have been obtained for 
the Burgers equation using an "instanton" approximationE3. 

Independently of the study of the stochastic models with a direction a considerable 
attention was paid in the last decade to investigation of models of random real asym- 
metric and complex non-Hermitian matrices. Eigenvalues of such matrices are generally 
speaking complex and so these models are quite different from models of random real sym- 
metric or Hermitian matrices. Starting from the first work in this directionil a number 
of publication contain discussion of properties of these models. Complex random 

matrices appear in study of dissipative quantum maps00 while real asymmetric random 
matrices have found applications in neural network dynamicsll@. Many interesting aspects 
of non-Hermitian matrices were discussed in preprintsil'S. Verv recently a new regime of a 
weak non-Hermiticity was found for complex random matrices^! In this regime an explicit 
formula for the density of complex eigenvalues was obtained by mapping the problem onto 
a zero-dimensional supermatrix a-model. 

Although one may guess that the models with the non-Hermitian or real asymmetric 
matrices should be related to disordered systems with non-Hermitian Hamiltonians, no con- 
vincing arguments have been given as yet. In fact, generally this is not true because, e.g. the 
models of open quantum dots described by non-Hermitian Hamiltonians can hardly corre- 
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spond to the models of random non-Hermitian matrices discussed in the literatureSIHI. 
However, as will be shown later, such a correspondence does exist in some limiting cases for 
the disorder models with a direction. 

The goal of the present publication is to develop a method that would allow to make 
analytical calculations for the disordered problems with a direction. This goal is achieved by 
modifying the supersymmetry technique in a way to include in the non-linear supermatrix 
ex-model terms corresponding to the imaginary vector potential. Although a proper cx-model 
for the physical real vector potential has been derived long agol, changing to the imaginary 
one is far from trivial and, as a result, a completely new term in the a-model appears. The 
zero-dimensional version of the a-model turns out to be exactly the same the one obtained 
in Ref.il for the model of weakly non-Hermitian random matrices. 

The supermatrix a-model derived below is valid in any dimension and can be a proper tool 
for studying the localization-delocalization transitions in one and two dimensions proposed in 
Reffi However, although one can use standard computational schemesi, the presence of new 
terms in the a-model make calculations with the known parametrizations of the supermatrix 
Q more difficult. Therefore, a new parametrization is suggested and corresponding Jacobians 
are calculated. To avoid "overloading" only zero-dimensional case is considered in this 
article. For the unitary ensemble the result of Ref.0 for the density of complex eigenvalues 
of weakly non-Hermitian random matrices is reproduced. The density function is a smooth 
function of the imaginary part of the eigenvalues, which shows that the probability of real 
eigenvalues is zero. 

In contrast, the density function for the orthogonal ensemble obtained below contains a 5- 
function, which shows that the fraction of states with real eigenvalues is finite. This is a new 
very unusual and interesting result. The entire function of the density of complex eigenvalues 
is obtained for the first time. In the limit of strong non-Hermiticity the probability functions 
for the both unitary and orthogonal ensembles correspond to the "elliptic law"@fi 

The basic results of this article have been presented in a short form elsewhereil. The 
article is organized as follows: 

In Section II models of disorder with a direction are introduced and their basic properties 
are discussed. Section III contains derivation of a supermatrix cr-model. In Section IV a joint 
probability density of complex eigenvalues is calculated for systems in a limited volume with 
broken time reversal symmetry (unitary ensemble). This is done by calculation of integrals 
over supermatrix Q for the unitary ensemble. A new parametrization for the supermatrices Q 
is introduced. In Section V similar calculations are carried out for the orthogonal ensemble. 
The result for the density of complex eigenvalues proves to be qualitatively different from 
that for the unitary ensemble. Section VI contains a discussion of the results obtained and 
comparison with some other works. In Appendix the Jacobians corresponding to the new 
parametrizations for the supermatrix Q are derived. 



II. THE MODEL AND ITS BASIC PROPERTIES 

The initial classical model of vortices in a (d + l)-dimensional superconductor with line 
defects considered in Refs.0@ contains an interaction between the vortices. In the corre- 
sponding quantum model of ^-dimensional boson this describes an interaction between the 
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bosons. The interaction is, in principle, very important. Its short range part does not allow 
bosons to condense at one localized state. At the same time if it is strong enough there 
can be only one boson in a localized state and the problem maps onto the model of non- 
interacting fermions. Of course, this is not true for extended states for which one should use 
the model of interacting bosons. 

It is clear that one should first understand which one particle states are localized and 
which are not. Therefore, as in Refs.@ii, it is reasonable to start with a rf-dimensional 
Hamiltonian H of non-interacting particles including a constant imaginary vector potential 
ih and random potential of impurities U (r) 

H = H + U(r), H = { V±i^l (2.1) 

2m 

where p = —iV and m is the mass of a particle (boson or fermion). 

The random potential U (r) is assumed to be distributed according to the Gaussian 
^-correlated law 

(U (r)> = 0, (U (r) U (r')> = -±-6 (r - r') (2.2) 

where r is the mean free time, v is the density of states at the energy e involved. As has 
been mentioned in the Introduction, the potential U (r) corresponds to the potential of the 
line defects and h to the component of the magnetic field for the model of the vortices. 
At the same time, the Hamiltonian H, Eq. ( |2.1| ) can describe other systems as well. So, 
we may study properties of the Hamiltonian H without recalling each time where it comes 
from. Some of possible applications of Eq. fl2.ip have been listed in the Introduction. The 
directed quantum hopping appears a new interesting possibility. The Hamiltonian Hl of a 
lattice version of Eq. (|2.1|) can be written as follows 



H L = ~ E E ( ehe X + +e,C r + e-^C+Cr+eJ +E f/ W (2.3) 
Z r v=\ r 

where c + and c are creation and annihilation operators and {e u } are the unit lattice vectors. 

Although Eq. ( |2.3| ) was used in Ref.0 only for numerical calculations, it has a clear 
physical application. It describes quantum hopping of a particle from site to site in the 
presence of a random potential. However, the hopping probability along h is higher than 
in the opposite direction. In other words, the Hamiltonian Hl describes a directed hopping 
in a random potential. The systems with the Hamiltonians H, Hl, Eqs. (|2.1|J2.3j ), are not 
invariant with respect to inversion of the coordinates even after averaging over impurities. 
At the same time, they are time reversal invariant and therefore essentially different from 
systems with real magnetic fields. 

If necessary the Hamiltonians H and Hl can be generalized to include the vector potential 
A corresponding to a physical magnetic field. This can be done by the standard replacement 

p - p--A (2.4) 

c 

in Eq. (|2.1|). Proper changes can also be done in Eq. 
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Of course, the vortex model of Refill corresponds to Eq. ( |2.1|) with A =0 but already 
the hopping model can be considered in an arbitrary magnetic field. Changing the magnetic 
field (or, more precisely, the vector potential A) results in a crossover between ensembles 
with different symmetries. In analogy with "conventional" (non-directed) disordered systems 
these ensembles will be called orthogonal and unitary. 

Although the Hamiltonians H and Hl, Eqs. (|2.1|, |2.3|) are not Hermitian this fact does not 
contradict to fundamental laws of nature. In the problem of the vortices in superconductors 
these Hamiltonians appear after a reduction of a (d + l)-dimensional classical problem to a 
c?-dimensional quantum one using the transfer matrix technique, which is a formal trick. As 
concerns the directed hopping model the vector h can appear as a result of a coupling with 
another system (reservoir) which is not necessarily in equilibrium. The latter system can be 
subjected e.g. to an electric field, there can be non-decaying currents in it, etc. Integrating 
out degrees of freedom related to the reservoir one obtains an effective Hamiltonian that 
does not need to be Hermitian. 

In other words, the non-Hermitian Hamiltonians appear at intermediate steps of calcu- 
lations and manipulations with them should be considered merely as formal computational 
tricks. The corresponding wave functions and eigenenergies are only formal objects as well. 
Of course, one should understand how to relate initial physical observables to quantities 
calculated with the non-Hermitian Hamiltonians. 

It is relevant to mention that a classical directed model that can be considered as the 
counterpart of the directed quantum problem has been introduced long ago§. This is the 
model of a directed percolation that can describe, e.g. spreading of infection or fire in 
a forest affected by wind. According to a discussion of Ref." critical behavior near the 
percolation transition in the model of the directed percolation is different from that of an 
isotropic model. The analysis of Ref.0 was based on a diagrammatic expansion. The bare 
Green functions (p) used in the expansion had the form 

G (0) (P) = 2 1 - (2.5) 

with a constant vector a. Comparing Eq. ( |2.5| ) with Eq. (|2.1| ) we see that is the Green 
function of the Hamiltonian Hq, which demonstrates that both models are really closely 
related to each other. 

Now, let us discuss following Refill basic properties of eigenstates of the Hamiltonian H, 
Eq. Q2.ip . Due to the non-Hermiticity of the Hamiltonian one should distinguish between 
right <pk (r) and left 4>k (r) eigenfunctions. They obey the following equations 

H<j) k (r) = e fc fc (r) , H T ^ k (r) = e fe fc (r) (2.6) 

where H T is obtained by transposition of the Hamiltonian H. For spinless particles the 
operation of the transposition means simply changing of the sign of the space derivative. 
The functions fa (r) are also called conjugate to 0& (r); for each eigenfunction one can 
construct its conjugate. The scalar product (fa, fai) of two eigenfunctions fa (r) and fa (r) 
is introduced as 

' h (r) fa, (r) dv (2.7) 
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Using Eq. (|2.7|) one can prove in a standard way the orthogonality of eigenf unctions corre- 
sponding to different eigenenergies. Together with the normalization condition this can be 
written as 

<) k (r) fa (r) dr =5 kk > (2.8) 

The eigenenergy e k in both Eqs. (|2.6|) is the same. Eq. (|2.8| ) enables us to reproduce basic 
properties of conventional (Hermitian) quantum mechanics replacing everywhere complex 
conjugates 0* (r) of the functions cj> k (r) by the conjugates (p k (r). However, the eigenenergies 
6 k in the non-Hermitian quantum mechanics are not necessarily real. They must be real only 
if the functions <p* k (r) and <j) k (r) coincide. In order to obtain well defined wave functions in 
the thermodynamic limit it is convenient to impose periodic boundary conditions. 

To understand better how the wave functions look like in different situations it is in- 
structive to consider a localized state with a localization center at a point Xq and extended 
states in the absence of impurities (for simplicity we may restrict ourselves with the purely 
one dimensional case). Assume that for h = the eigenf unctions (f> k and the eigenvalues 
ej^ are known. Then, the functions 

fc (x) = e hx <$ ] (x) , 0, = e-^f (x) (2.9) 



(0) 
k ■ 



are solutions of Eqs. fl2.6|) with the eigenenergy e 

At the same time, in order to satisfy the boundary conditions the function <p k and 4> k 
may not grow. If the function <f> k ^ (x) is exponentially localized at a distance l c , the function 
(j) k (x) takes the form 

4>k (x) = C exp (h (x — Xq) — l^ 1 \x — x \j (2-10) 



The function <p k (x), Eq. (|2.10|) , and the corresponding function (p k (x) does not grow at \x\ — > 
oo only if \h\ < The point \h\ = Z" 1 was identifiedEl with a localization-delocalization 
transition. 



In the region \h\ > l c 1 the functions <p k given by Eqs. (|2.9| , p . 1 0|) are longer eigenfunctions 



because they do not satisfy the boundary conditions. To get an idea how the eigenfunctions 
look like in this region we may neglect the disorder potential. Then, the plane waves 

fc = L~ 1/2 e ikx , fc = L~ 1/2 e~ ikx (2.11) 



where L is the length of the sample, are proper solutions of Eqs. ( |2.6j ) satisfying the boundary 
conditions. However, in this case the eigenvalue e k is no longer real 

* = ^2 (2.12) 
2m 

We see that the question about whether an eigenfunction in the presence of the imaginary 
vector potential is localized or extended is closely related in the thermodynamic limit to 
the question whether the corresponding eigenenergy is real or complex. The arguments 
presented are qualitative but they were confirmed by numerical calculations^. 
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It is clear from the previous discussion that it is very important to understand when 
eigenenergies are real and when they become complex. A convenient function characterizing 
the system is the joint probability density of complex eigenenergies P (e, y) defined as 

P(e:V) = ^(j:He-e' k )5(y-e'Q\ (2.13) 

where e' k and e k are the real and imaginary parts of the eigenenergy e k , V is the volume and 
the angle brackets stand for averaging over impurities. If all states are localized, such that 
e' k = 0, the function P (e, y) equals 

P(e,y) = u(e)5(y) (2.14) 

where v (e) is the average density of states. 

If all states are extended the function P (e, y) should be a smooth function of both 
variables. In some cases physical quantities can be expressed directly through the function 
P(e,y) although other correlation functions are also of interest. The rest of this article 
is devoted to reduction of the function P(e,y), which is the simplest non-trivial function 
characterizing the system, to a correlation function in a supersymmetric cx-model and to 
some calculations with this model. This is the first attempt of a quantitative analytical 
study of the disordered directed quantum systems. 



III. DERIVATION OF cr-MODEL 

According to the standard procedure of derivation of the supermatrix cx-mo defli one 
should express the physical quantity in terms of retarded Gf and advanced Gf Green func- 
tions of the Hamiltonian. Usually the average density of states that can be expressed through 
the average of one Green function is not an interesting quantity because it does not distin- 
guish between localized and extended states. The density of complex eigenvalues P (e, y) 
is definitely more interesting but how to express it in terms of integrals over supervectors, 
which is the first step of derivation of the cr-model? 

The problem is that it is not clear how to write the function P (e, y) in terms of the 
functions Gf, Gf. However, even if this representation existed it would not help. Using the 
spectral expansion of the functions G R,A 

G ^ {r y ) = E M^Ml (3.1) 
^ e - e k ±ib 

we see that if some eigenenergies e k are complex the function Gf (Gf^j is no longer analytical 
in the upper (lower) half plane of complex e. But the very possibility to rewrite the Green 
functions in terms of convergent Gaussian integrals over the supervectors was based on the 
assumption that the eigenenergies were real. 
Another possibility is based on the relation 

6 (a) 6{b) = l Km - f (3.2) 
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that holds for real a and b. With Eq. (|3.2|) the density function P (e, y) can be rewritten as 

P ^ V) =^F Km (E [( e - + (V ~ 4f + 7 2 ] " 2 ) (3.3) 
Using the orthogonality of the eigenfunctions 4> k Eq. (|3.3| ) can be also represented as 

P ( e = lim / 5 (r, r') 5 (r\ r) drrfr' (3.4) 
where the function B (r, r') has the form 

B(r,r') = E 7 t ^ (r0 2 (3.5) 

The representation of the density function P (e, y) in by Eq. ( |3.4|) is very convenient because 
it allows to rewrite this function in terms of a Gaussian integral over supervectors. 
In order to derive a proper expression let us introduce an Hermitian operator M 

M=( .f™7 £ , HH "r V \) (3-6) 



-i (H" - y) - (H 1 - e) 



where 



H>=±(H + H+), H" = - l -(H-H + ) (3.7) 

In Eq. (|3.7| ) the symbol "+" means Hermitian conjugation. For real Hamiltonians this 
conjugation coincides with the transposition "T" . However, let us write formulae in a general 
form such that the Hamiltonian H may include magnetic interactions and be complex. 

Instead of manipulating with the non-Hermitian operator H one can try to use the 
Hermitian operator M. To follow the standard procedure of the supersymmetry technique 
one should find first the eigenstates of this operator. For the complex non-Hermitian operator 
H one can write 4 equations for the eigenstates 

H<p k = e k 4> k , H T <j) k = e k <f) k (3.8) 



H*ct>l = e* k </>i, H+ti = e\4>l (3.9) 



Eqs. ( |3.9| ) are merely complex conjugates of Eqs. ( |3.8| ). 

Now, let us introduce two sets of 2-component vectors u k and v k 

u k = I ( 4> k + <P* k 4>k - <t>l ) ; v k = - ( (j) k - <p* k (j) k + <p* k ) 
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Using the orthogonality of the eigenfunctions fa, Eq. fl2.8|) , one can prove the orthogonality 
of the vectors u k and v k 

j u k (r) uy (r) dr = J v k (r) Vy (r) dr =5 k y (3.11) 

J u k (r) v k > (r) dr = J v k (r) uy (r) dr =0 

It is not difficult to see that the vectors u k (r) and v k (r) are eigenvectors of the matrix 
operator M satisfying the equations 

Mu k = M k u k , Mv k = M k v k (3.12) 

where the matrix M k equals 

*-(4rv ( &~4) (3 - l3) 

and e' k , e'J, are the real and imaginary parts of the eigenenergies e k . 
Using the identity 

\Tr (M k + Z7)" 1 = -y— L -y— - (3.14) 

one can see that the functions B(r, r'), Eq. (|3.5|) , are closely related to the operator M. 
The only thing that remains to be done is to express the matrix (M k + and then the 

operator (jvl + i^j in terms of a Gaussian integral over supervectors. 

The operator M is Hermititian, its eigenvectors u k and v k , Eqs. (|3.10|) are known and 



therefore we can follow the standard procedure of the derivationQfl. Changing from the 
Hamiltonian H to the operator M we had to double the size of the relevant matrices. 
This means that in order to write proper Gaussian integrals we should use, as usually, 8- 
component supervectors ip (r). In fact, one comes to supervectors if) with exactly the same 
structure as previously^! 

/ „om \ 1 / m* \ 1 / am* \ 

V» m = !m , # m = ~m Im > = -= % m ) (3.15) 



m = 1,2; x" 1 an d 5" m are anticommuting and commuting variables respectively. 

Let us present several important intermediate steps of the reduction of the operator 

(M + , to the functional integral over if) (r). First, we have 

(ry + MkY 1 = -i J (a k a* k + b k b* k ) exp (-L k ) dR k = -i J (a k a* k + p k p* k ) exp (-L k ) dR k 

(3.16) 

where a k , b k and a k , p k are commuting and anticommuting variables, respectively, dR k stands 
for the elementary volume in the space of these variables. The function L k in Eq. (|3.16|) 
equals 
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L k = -i(a% bl)(i 1 + M h )^-i(a* k pi) (vy + M k ) ^ (3.17) 



The vector fields x ( r ) an d S (r) are introduced as 



*(')=[ *2 (r) ) = E i^u k (r) + b k v k (r)) , (3.18) 



^0) = ( £2 (r) ^ = E (^fe ( r ) + ( r )) 



where the vectors and are defined in Eqs. (|3.10|) . 

With these definitions one can express the functions B(r,r'), Eq. (|3.5| ), in terms of 
Gaussian integrals over the vector fields x(r) an d S(r). The derivation is based on the 
identity 

J S* (r) MS(r) dv=Y^{alb%)M k ^ £ j (3.19) 

that can be proven using Eqs. ( |2.8|J3.10| ) (the same for x(r)). Less trivial is the expression 

ij [ S* (r) S (r) dr (3.20) 



Using the expansion, Eq. (|3.18|) , we can see that the integral, Eq. (|3.20|) , contains non- 
diagonal with respect to k, k' terms. For example, there is the following term 

^J» fc / ^(r)0 fc (r)dr (3.21) 

Z k,k' J 

For Hermitian Hamiltonians the integral in Eq. ( |3.21| ) would give 5 k y. However, generally 
it is not zero for arbitrary k and k' because the orthogonality relation, Eq. (|2.8|) , contains 
<p k but not <p* k . Fortunately, this does not create difficulties in the limit of small "vector 
potential" h that is of the main interest in the present work, because the difference between 
(fi k and 4>* k is small. This allows us to write 

5> fc = f\r (r) (i 7 + M) x (r) + S* (r) (i 7 + M) S \v)]dv (3.22) 
k J 

Although one can use Eq. (|3.22|) as an effective Lagrangian, it is convenientSi to unify all 
components of the vectors x, x*, S, and S* into the supervector ip of the form, Eq. ( |3. 15[ ) . 
As a result, one comes to integration with the weight exp (— C), where the Lagrangian £ 
takes the form 

C = -i[i)(r)(Ha + U(r))i)(r)dr (3.23) 

where the "charge-conjugate" supervector ip (r) is the same as in Refs.ii. The 8x8 matrix 
operator Ho can be written as 
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'Ho — T^oo + T~toi > 



(3.24) 



n 



oo 



H' -e + ijA, H 01 = iA 1 {H" + yT 3 ) 



In the continuum model, the "imaginary" part H" of the Hamiltonian H, Eqs. (|2.1] , |3.7D 
has the form 



m 



(3.25) 



The diagonal matrices A and r 3 are the same as in Refs.Qfl The matrix Ai anticommutes 
with the matrix A and also consists of unit 4x4 blocks. The explicit form of these matrices 
is 



A 




Ai 




^3.26) 



Eq. ( |3.23|) is similar to the corresponding equation for localization problemsLrO and in the 
absence of 7ioi these equations would coincide. All new physics comes from the operator 
Hoi- A magnetic field can be included into 7i o i n a standard way. 

All subsequent manipulations are the same as in Refs.0!. First, one averages over the ran- 
dom potential U (r) using Eq. (|2.2|) and comes instead of Eq. ( |3.23| ) to a regular Lagrangian 
C 



(r) (r) + 



Ahvt 



(if> (r) i/j i 



dr 



(3.27) 



Then, one decouples the interaction term in Eq. ( |3.27[ ) by integration over a supermatrix 
Q and integrates over the supervector ip assuming that the supermatrix Q varies in space 
slowly. After that one comes to an integral over Q with the weight exp (-F [Q]). The 
functional integral over Q is calculated using a saddle-point approximation. At the saddle- 
point the supermatrix Q does not depend on coordinates and in the limit of small 7ioi and 
7 one obtains the standard equation 



Q(r) 



1 

iris 



-iUno + 



Q(r) 
2r 



(3.28) 



which leads to the constraint Q 2 = 1. Now, one has to expand the free energy functional 
F [Q] near the saddle-point in iJ i> 7 an d VQ. As a result the functional F [Q] acquires the 
form of a cr-model 

f iq] = yJ STr[D ° (vg + h [g ' Ai])2 " 4 (7A + yAlTi) Q]dr (3 - 29) 

where D is the classical diffusion coefficient, [.,.] is commutator and STr stands for su- 
pertrace. Eq. (|3.29| ) is written in the absence of a magnetic field. The expansion near 
the saddle-point leading to Eq. ( |3.29| ) is justified provided y <^ r _1 and h l" 1 , where 
I is the mean free path. The supermatrices Q are the same as those for the orthogonal 
ensemble&i. This case can correspond to the problem of vortices in superconductors with 



12 



line defectslll. If for some other problems one has to include in the Hamiltonians H, Hl 
the physical vector-potential A corresponding to a magnetic field, the standard derivation 
shows that the proper cx-model is obtained from Eq. (|3.29|) by the replacement 



26 

V<2^ VQ A [Q, ra] (3.30) 

c 

In the limit of a strong magnetic field one can neglect fluctuations of a certain symme- 
try (cooperons). Then, Eq. Q3.29| ) is still valid but the supermatrices Q should have the 
symmetry corresponding to the unitary ensemble. 

The free energy functional F [Q], Eq. ( |j.29j ) has two additional with respect to the 
functional used for "conventional" disorder problems. These terms contain the matrix A 1; 
which leads to new effective "external fields" in the free energy. We see from Eqs. (|3.29|j3~30| ) 



that h and A enter F [Q] in a different way. A simple replacement A -^ih in the cx-model 
of Refs. (i'i) would give a wrong result. This reflects the fact that a no n- zero A violates 
the time reversal symmetry while h can break only the symmetry with respect to inversion 
of coordinates. 

In order to write express the density function P(e,y), Eq. ( |3.4j , |3.5| ), in terms of a 
functional integral over Q one should know not only the weight exp (-F [Q]) but also a 
pre-exponential functional A [Q] . It can be derived from Eqs. ( |3~4|,|3.5|) in a standard way. 
One of the functions B can be written using the first line of Eq. (|3.16|) and the other using 
the second one. As a result, one obtains in the pre-exponential a product of four different 
components of the supervector if>; two of them are at the point r while the other two at the 
point r'. After averaging over the random potential U (r) and decoupling of the effective 
interaction in Eq. ( |3.27| ) by integration over the supermatrix Q one has to compute Gaussian 
integrals over This can be done using the Wick theorem. In the limit r _1 <C {yV) 1 one 
may take into account only pairing of two if) at coinciding points. The rest of the calculation 
is simple and one obtains 

P(e, 2 /) = -lim^ r [A[Q]exp(-F[Q])dQ, (3.31) 

where 

A [Q] = J [(Qll (r) + QH (r)) (Qll (r') + Qf A (r')) (3.32) 

- (Qll (r) + Qll (r)) (Qll (r') + Q% (r')))drdv> 

Numeration of the matrix elements in Eq. ( |3.32| ) is standard!!!. 

Eqs. (|3.29| - |3~32"|) solve the problem of mapping of the density of complex eigenvalues for 



disorder models with a direction onto a supermatrix cx-model. The density function P (e, y) 
depends on the real part e of the eigenenergies through the parameters v and D that are 
dependent on e. The dependence on the imaginary part y is more complicated. Remark- 
ably, the cx-model derived differs from the cx-model for localization problems by additional 
"external fields" only. This simplifies calculations because one can use well developed com- 
putational schemes. 
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The cr-model, Eqs. ( p.29[ - p73^ ) can be used in any dimension. The one- dimensional 



version describes "quantum wires" or, in the language of the superconductor model, to 
vortices in a slab. According to a discussion of Ref Jia, in one-dimensional models there has to 
be a localization-delocalization transition. If this is true for thick wires the one-dimensional 
cr-model should undergo a phase transition when changing the value of h. However, study 
of the one-dimensional model is more difficult than of the zero-dimensional one. Leaving 
higher dimensional problems for future investigation let us concentrate in the next Section 
on calculating the density function P (e, y) for a sample with a finite volume. This situation 
is described by the zero-dimensional cr-model. 

IV. DENSITY OF COMPLEX EIGENVALUES IN A LIMITED VOLUME: 

UNITARY ENSEMBLE 

If disorder is not very strong there is a regime when physical quantities can be obtained 
from the zero-dimensional {OD) cr-model. This is the limiting case when one considers only 
supermatrices Q that do not vary in space. For the problem of level statistics in Hermitian 
models the OD cr-model is obtained in the limit uo <C E c , where E c = tt 2 Dq/L 2 is the 
Thouless energy {L is the sample size)SH. If the sample is connected with leads and the 
energy levels are smeared the OD case is possible provided the level width does not exceed 
E c . If the disorder is strong or the sample has one- or two-dimensional geometry, such that 
the localization length L c is smaller than the sample size, the OD limit cannot be achieved. 

It is clear that the situation with the directed problems involved should be similar and 
one can come to the OD cr-model provided h, y, and 7 in Eq. ( |3.29| ) are not very large and 
disorder is not very strong. For the model of vortices in a superconductor the OD limit for 
the cr-model would correspond to a sample with a finite cross-section perpendicular to the 
line defects. 

Neglecting all non-zero space harmonics in the free energy functional F [Q] one can 
rewrite Eq. ( p. 29 ) as follows 



F [Q] = STr [Q, A,] 2 - ^A^Q - 1aq\ (4.1) 

where 



and A=(uV) 1 is the mean level spacing. 

The distribution function P (e,y), Eqs. ( j3.31] , j3.32"D , takes the form 



P (e, V) = I™, / A [Q] exp (-F [Q]) dQ, (4.3) 

A [Q] — (qII + Q42) + QlT) ~ (Q42 + Q42) (Q24 + Q 

with F [Q] determined by Eq. (|4.1|) . 
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The non-zero space harmonics can be neglected provided the following inequalities are 
fulfilled 



y < E c , 7 < E c , h<^L 1 (4.4) 

where L is the sample size. 

To obtain the function P (e, y) one should calculate in Eq. ( ^4.3| ) a definite integral over 
the supermatrices Q. The structure of supermatrices Q is the same as in Refs.Qfl and, in 
principle, the way how to compute the integral is clear. As usual, all manipulations are 
simpler for the unitary ensemble and therefore let us start with this case. 

However, already before an explicit calculation of the integral in Eq. Q4.3| ) an interesting 
observation can be made. We know that the OD version of the a-model for Hermitian 
disordered systems can also be derived from random matrix mo dels!. In fact, it is the way 
how the equivalence of between disordered systems in a limited volume and random matrix 
theory (RMT) was finally established. Now, a natural question arises: do the random models 
with a direction considered in the present work correspond to a RMT? 

Of course, this cannot be a model of Hermitian or real symmetric matrices because 
in this case all eigenvalues must be real. So, one should think of ensembles of random 
real asymmetric or complex non-Hermitian matrices. Study of random complex matrices 
without the requirement of Hermiticity has started quite long agoil and since then models 
of non-Hermitian or real asymmetric random matrices have been considered in a number 
of publicationsiB@0iHI. The ensembles of real symmetric random matrices have found 
applications in e.g. neural network dynamicsil'll while the ensembles of complex random 
matrices appear in study of dissipative quantum mapsHll. One of results obtained is that, 
for Gaussian ensembles in the limit of a large size N of the matrices, the eigenvalues are 
uniformly distributed in an ellipseiSi. 

Recently, an ensemble of "weakly non-Hermitian" random matrices X was introduced^. 
It was assumed that these matrices had the form 

X = A + iaN- 1/2 B (4.5) 

with N x N statistically independent Hermitian matrices A and B, and a number a of the 
order of unity. The matrices A and B obeyed Gaussian distributions with the probability 
densities 

<x exp (- ^rri 2 ) , V (6) oc exp (- AjyB 2 ) (4.6) 



V(A) 



where J has order of unity. 

The parameter aiV~ 1//2 is a measure of the non- Hermiticity and is always small for N — > 
oo and a finite. The authors of Ref.0 calculated a density of complex eigenvalues similar to 
the function P (e, y), Eq. ( |2.13| ) and demonstrated that this function has a finite limit when 
iV —>■ oo. At the same time they did not point out any direct physical applications. For 
computation of the function P (e, y) they used the supersymmetry technique. Remarkably, 
a cx-model derived in Ref.0 is exactly the same (although numeration of elements of the 
matrix Q is somewhat different) as the unitary version of OD a-model, Eq. ( f4.1| ). The pre- 
exponential is different but this is natural because another (less direct) way of calculating 
the function P (e, y) was used. 
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The same form of the cx-model obtained for these two different models shows that the di- 
rected disordered model with broken time-reversal invariance in a finite volume is equivalent 
to the model of weakly non-Hermitian matrices. Apparently, the same equivalence holds be- 
tween the time reversal invariant model of disorder and models of weakly non-symmetric real 
matrices. However, it is relevant to emphasize that not every non-Hermitian Hamiltonian 
corresponds to the models of non-Hermitian or non-symmetric real matrices. For example, 
models of open chaotic billiards are described by Hamiltonians with additional imaginary 
terms (see, e.gJji). These Hamiltonians do not seem to be equivalent to the random matrix 
models of RefM 

Now let us show how explicit calculations in Eqs. ( |4.1| , [4.3[ ) can be performed. First 



of all one should choose a proper parametrization of the supermatrices Q. The authors 
of Ref.0 used the parametrization of Ref.0 ("standard parametrization" in terminology of 
Ref.i). This parametrization has been used for solving many interesting problems. However, 
due to presence of the new terms in the free energy F [Q], Eq. ( |4.1| ), this parametrization is 
not as convenient as before! because now F [Q] would contain not only the "eigenvalues" 6 
but also many other variables. 

As concerns the unitary ensemble, the computation of the function P (e, y)is still possible 
although is very lengthy^. At the same time, calculations for the orthogonal case using the 
standard parametrization do not seem to be possible at all due to unsurmountable technical 
problems. 

Fortunately, one more parametrization is possible that is perfectly suitable for the present 
problem. To some extent it resembles the parametrization used to study the crossover 
between the orthogonal and unitary ensemblesEafl. Of course, it should be written for the 
orthogonal and unitary ensembles in a different way but the main structure is the same. 
Let us show in this Section how the function P (e, y) can obtained for the unitary ensemble 
using this new parametrization (It can be named "non-Hermitian parametrization"). The 
orthogonal ensemble will be considered in the next Section. 

The supermatrix Q in the non-Hermitian parametrization is written in the form 

Q = TQqT (4.7) 

where T should be chosen to satisfy the relations [T, Ai] = 0, TT = 1. The bar stands for 
the "charge conjugation" defined in RefsBi. It is clear that with such a choice the function 
F [Q] would depend on Q only (for the unitary ensemble one has also [Q , r 3 ] = 0). 
The central part Qq in Eq. (|4.7|) is taken in the form 



-T 3 Sill ip — COS (fi 

while the supermatrix T can be chosen as 





r -Ml-^ (4 - 9) 

The supermatrices 9, u, v are equal to 
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_ / 1 - 2r]f} 2r] \ _ f 1 - 2kk 2k \ 

U ~ \ -2fj 1 - 2^77 J ' V ~ \ -2k \-2kk) 

The 2x2 matrices p, x, 9, an d Q\ are proportional to the unit matrix, the matrices rj, k are 

-(;'_:.)■ -(0-°.) ^ 

where 77, 77*, k, and k* are anticommuting variables. The conjugate matrices fj and k are 
the same as in RefsBi. To understand better the structure of the supermatrix Q given by 
Eqs. ( |4.7| - |4.11]) it is instructive to write it neglecting all Grassmann variables. Then, one 
can write separately the compact and noncompact sectors. The compact sector takes the 
form 



cos 9 cos (p — T3 simp + i sin 9 cos <p 

-T3 simp — i sin9 cosp> — cos 9 cosp 



(4.12) 



whereas the noncompact sector is written as 

(4.13) 



cosh 6*i cos x —i T 3 sinh x — sinh 9\ cosh x 

-itz sinh x + sinh 9\ cosh x ~ cosh. 9\ cos x 



Comparing Eqs. ( |4.12| , |4.13D with the corresponding expressions for the supermatrix Q in 
the standard parametrizationl'i one can understand that in order to specify the supermatrix 
Q unambiguously the following inequalities should be imposed 

- 00 < % < 00, -00 < 9i < 00, -7r < 9 < 7r, — 7r/2 < p < 7r/2 (4.14) 

To start computation with the parametrization, Eqs. ( [4.7| - f4.11] ), one should derive first the 
proper Jacobian. The derivation is presented in the Appendix. The final result for the 
elementary volume [dQ] reads 

[dQ] = J^JgdRsdRF, dRs = d9d9\dpdx, dRp = d^drfdndn (4-15) 

where 

_ 1 cosy cosh x f4 s 
o7T [sinh x + 1 sm V 9 ) 



32tt sinh 2 \ {9 X + it 



(4.17) 



Substututing Eqs. ( |4.7H4.11| ) for Q in Eq. ( |4.1|) one can rewrite the function F [Q] in the 
limit 7 — ► as 
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F [Q] = a 2 (sinh 2 x + sin 2 ip 



ix (sinh \ + i sin if) 



(4.18) 



(The limit 7 — > is taken in the beginning of the calculations because in the present 
parametrization this does not lead to additional convergence problems). The function F [Q], 
Eq. ( fl.l8j ), does not contain the anticommuting variables and therefore one can easily 
integrate over the supermatrix u. Writing in Eq. (|4.3|) the supermatrix Q, Eqs. Q4.7| - [4.9|) , 

as 



Q = uQu 

with u from Eq. ( |4 . 1 0| ) and integrating over rj, rf one obtains 



(4.19) 



Plr.u) = ^J (5Tr(r 3 A 1 g)) 2 exp(-F[g])rfg 



A.-KV d 2 
A dx 2 



exp(-F [Q])dQ 



(4.20) 



where the elementary volume dQ differs from [dQ] by the replacement dRp — > dRp = 
dtzdtz* and F[Q] is given by Eq. (|4.18|) . Although Eq. ( |4.20| ) is quite simple, one more 
difficulty should be overcome. The problem is that the integrand in Eq. fl4.20|) does not 
contain the variables k, k* and, at first glance, the integral must turn to zero. However, the 
Jacobian J e , Eq. ( f4.17[) , is singular for 9, 61 — > and this singularity is not compensated by 
the integrand. So, one obtains an expression of the type x 00, which is a usual phenomenon. 
Different procedures how to make the integral well defined have been worked out (for a 
detailed discussion see!). The simplest way is to rewrite Eq. (|4.3|) as 



where 



P (e, y) = P m (e, y) - J A [Q] (exp (-F [Q]) - exp (-F m [Q])) dQ (4.21) 



P m (e, y) = -^J A [Q] exp {-F m [Q]) dQ, 



(4.22) 



F m [Q] -F[Q} = -mSTr (TATA) = -mSTr (fiA) 
The supermatrix Q in Eq. (|4.22|) can be chosen as 

tt = TAT, T = uf 

The parameter c in Eqs. Q4.21 , 4.22j) is arbitrary. Using Eq. ( [4.9|) we see that 

- STr (TAT A) = 4 (cosh#i - cos^) 



and thus, the singularity at d\ = 9 = coming from the Jacobian in Eq. ( [4.21 



(4.23) 



(4.24) 



is com- 



pensated by the integrand. After integration over 77, 77* the integrand does not contain the 
anticommuting variables k, k* and the integral vanishes. Therefore, the function P m (e,y), 
Eq. Q4.22D , does not depend on m and one can calculate the integral in the limit m — > 00. 
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In this limit only small deviations of the supermatrix from A are essential. Using the 
representation, 

n = A{l + iW){l-iW)-\ W=(yl B Q ),B=( a _ "A (4.25) 

expanding O in W up to quadratic terms and calculating the Jacobian in this approximation 
one can see that in the limit m — > oo 

J exp (-mSTr (OA)) rfO = 1 (4.26) 

The supermatrix T can also be represented through W and calculating the corresponding 
Jacobian one may expand up to quadratic in W terms. As concerns Q in the other terms 
in the integrand in Eq. ( |4.22| ), one should replace in the limit m —>■ oo the supermatrices T 
by 1. One can check also that now the Jacobian of the transformation from the matrices T 
and u to T equals to —1 and not to Jq as it was with the initial parametrization for T, Eq. 

(ID- 

So, calculating the integral, Eq. (|4.3|) , one should replace the supermatrix T in the 
integrand by 1. In the elementary volume [dQ], Eq. fl4.15| ), one should omit the multiplier 
Jodndtz* and change the sign of the rest. 

As a result of all these manipulations one comes to the following expression for the 
function P (e, y) 

P (e, y) = ~^J ( STr (^AiQo)) 2 exp (-F [Q„]) J v dipd X (4.27) 



with Qq from Eq. (|4.8|) and J v from Eq. ( |4.16 ). The function F [Q ] is given by the R.H.S 



of Eq. ( |4.18| ). The limits of integration over tp and x are determined in Eqs. ( |4.14| ). 

The further calculation in Eq. ( 4.27|) is very simple because the function in the pre- 



exponential is proportional to J^ 1 . Changing the variables of integration z = sinhx, t 
sin (p, one is to calculate a Gaussian integral over z, and the final expression takes the form 

P (e, y) = V -^- exp i^^^j cosh xt exp (-a 2 t 2 ) dt (4.28) 

The function P (e, y) is properly normalized and one obtains using Eq. ( |4.2| ) 

P(e,y)dy = l (4.29) 



The density of complex eigenvalues P (e, y), Eq. (|4.28 ), agrees precisely with the correspond- 



ing function for weakly non-Hermitian random matrices obtained in Ref.0. The parameters 



a and A are related in this case to the parameters in Eq. (4.5 



as 



a = V2tvJu (e) a, A = (u (e) A^ 1 , v{e) = (2njy 1 \fi - (e/J) 2 (4.30) 
and x = 2tip (e) yN . 
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The agreement can serve as a proof of the equivalence between the directed disorder 
models in a finite volume (with broken time reversal invariance) and the models of non- 
Hermitian matrices defined by Eqs. (|4.5| , |4.6| ). The function P (e,y) is represented in Fig.[TJ. 



Its basic properties have been discussed in RefHO. 

The density of complex eigenvalues is a smooth function at any finite a, which means that 
any finite non-Hermiticity smears all eigenenergies making them complex. The probability 
of real eigenvalues is negligible. For a ^> 1 the integral in Eq. Q4.28p can be calculated 
analytically using the saddle-point method. In the interval \x\ < 2a 2 the integrand as a 
function of t has a sharp maximum in the domain of the integration and the integral can be 
extended to infinity. For \x\ > 2a 2 the function P decays fast. As a result one obtains 

F ^ V >~ 2a* A 0, \x\>2a 2 ^ dl) 



Eq. ( 4.31 ) shows that for a ^> 1 the density of imaginary parts y of eigenvalues at a fixed 



real part is homogeneous in the interval x G (—2a 2 , 2a 2 ). Using Eq. (|4.30|) for v (e) and a 
we can rewrite the result expressed by Eq. ( f4.31p in terms of distribution of eigenvalues in 
the complex plane. In such a formulation, Eq. (|4.31|) means that the complex eigenvalues 
are distributed homogeneously within the ellipse 



e ^ 2 



This is the "elliptic law" found in Refs. , which is natural because the limit o> 1 should 
correspond to a "strong" non-Hermiticity. At the same time, it is clear the elliptic law is 
model dependent. For the models of disorder considered in the present paper the density of 
complex states essentially depends on y only. 

In the opposite limit a«l the density of complex states P (e, y) takes the form 

V\Fti ( X 2 ^ 



P( £ , s ).-^ e xp(--j (4.33) 

The Gaussian form of the function P can be easily understood starting from the random 
matrix model, Eqs. ( |4.5| , [4.6|) . The function P (e,y) can be written as 

N 

P (e, y) = iV- 1 £ (5 (e - e'J 8 (y - <)) (4.34) 



n=l 

A' 

2tiN 



I N roo 

^ E dkelkV ( e " e n) exp (-ifc£)> 



Where the angle brackets (...) stand for the averaging over the matrices A and B, Eq. ( [4.6|) . 
In the limit of small a the imaginary part e" can be obtained using the standard perturbation 
theory. In the first order one has 

el = $* m B$ m (4.35) 
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where m is the eigenvector of the matrix A corresponding to the eigenvalue e' m . Substituting 
Eq. (|4.35|) into Eq. (|4.34j ) one can immediately average over the matrix B. Using the 
orthogonality of the eigenvectors <j) m one can write the result of the averaging as 



1 N 
— Y 



y°° dke iky Id (e - e'J exp 




(4.36) 



...) stands for averaging over A. Integrating over and using Eq. ( |4.30| ) one comes 

'dm 



where 

to Eq. ([4.33 ). As concerns the models of disorder, Eqs. (|2~T -2.3), even the asymptotics, 
Eqs. ( |4.31| , |4.36| ), have not been known before and it not clear how to reproduce them using 
simple arguments. 

Are the result obtained in this Section general and one cannot expect anything new for 
the orthogonal ensemble? Of course, there is no reason to hope that Eq. (|4.28| ) describes 
the orthogonal ensemble as well but are the asymptotics in the limits a ^> 1 and a < 1, 
Eqs. ( ^311 - ^331) still correct? 

The orthogonal ensemble of random matrices can be introduced again by Eqs. ( f4.5| , |4.6|) 
but now the matrices A and B should be real symmetric and antisymmetric, respectively. 
One should also make in Eq. (|4.5| ) the replacement a — * — ia. As concerns the asymptotics 
in the limit a ^> 1 the same elliptic law as in Eq. ( [4.32[ ) has been recoveredil. At the same 
time, one can expect completely different behavior for a < 1. This can be seen easily from 
the fact that the first order of the perturbation theory corresponding to Eq. ( |4.35| ) gives 
zero and one cannot derive Eq. ( |4.33[ ) as before. In fact, the density of complex eigenvalues 
P (e, y) is singular at y = 0. Study of the orthogonal ensemble is presented in the next 
Section. 



V. DENSITY OF COMPLEX EIGENVALUES IN A LIMITED VOLUME: 

ORTHOGONAL ENSEMBLE 

To compute the density of complex eigenvalues P (e, y) for the orthogonal ensemble one 
can start as previously from Eqs. ( [4. 1| - |4.4| ) but now one should use supermatrices Q with 
the structure corresponding to this case. As has been mentioned, the presence in Eq. ( |4.1| ) 
of the new term with the matrix Ai makes the calculation very difficult even for the unitary 
ensemble and hardly feasible at all for the orthogonal one. So, as in the preceding Section 
a new parametrization for Q should be designed. 

Let us write the supermatrix Q in the form 

Q = ZQ Z, Z = TY (5.1) 



with the supermatrices Qo and T specified by Eqs. (|4.8|44.1ll ) and choose the supermatrix 
Y as follows 

Y = Y RS, Y = Y 3 Y 2 Y 1 (5.2) 
The supermatrix Y\ entering Eq. (p.2|) is 
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w 
w 



w 



w 
1 



w 



cos(/i/2) — sin(/i/2)\ 
sin(/i/2) cos(yu/2) J 



(5.3) 



The supermatrix Y 2 is equal to 

Y 2 = 



cos (4/2) -i sin (4/2 
-i sin (§ 2 /2) cos (4/2 



5 ^2 



The supermatrix Y3 is 
r , exp(z/3/2 









i9 2 Ti 



/3r 3 
/3ir 3 



T 3 



1 

1 



1 
-1 



(5.4) 



exp (i/3/2) y 

The supermatrices i? and S contain remaining Grassmann variables and are written as 



(5.5) 



R 



R 
R 



R 



1 - 2pp 2p 
-2p 1 + 2pp 



P 



and 



S 



1 - 2<t 2 2ia 
2ia 1 - 2<j 2 



a 



(7 

a 



a 



P 
-p* 



a 

-(7* 



(5.6) 



(5.7) 



where p and a are conjugate to p and o. 

The parametrization for F, Eqs. ( |5.2| - |5.71 ), is chosen in such a way that \Y, Ai] = 0. To 
specify the supermatrix Q unambiguously one should restrict variations of the variables by 
certain intervals. This can be done as the preceding Section by comparing the bosonic "skele- 



ton" of Q written in the parametrization, Eq. ( |5.1[ - |5.7| ), (let us called it "non-symmetric 
parametrization" ) with the standard parametrization of Refs.Bi. As a result one can write 
the following inequalities 



< % < OO, -7T/2 < if < 7r/2, -OO < 9 1 < OO, -7T < 9 < 7T 
< 6 2 < OO, < p < 7T, < (3 < 7T, < f3i < 2-K 



(5. 



The next step is to calculate the Jacobian. The derivation is presented in the Appendix and 
the final result for the elementary volume [dQ] is 



[dQ] — J^JgJ ^JcdR^dRpdRiBdRip 



(5.9) 



In Eq. ( |5.9|) , J v , Je,dRB and dRp are given by Eqs. G4.15| - [4.17|) . The additional quantities 
entering Eq. ( |5.9| ) are equal to 



Ju 



sinh 8 2 sin p 



2 8 vr 2 (cosh6» 2 - cospY 



(5.10) 



Jr. 



4 sin 2 ip 



(sinhx — isinipY 



(5.11) 



22 



and 



dRiB = dpd9 2 dpdPi, dRip = dada* dpdp* (5-12) 
The free energy F [Q], Eq. ( |4.1| ), takes in the limit 7 — > the following form 

F [Q] = a 2 (sin 2 (p + sinh 2 x) + x [{cos p sin (p — i cosh 82 sinh x) (5.13) 

+4 (era* + pp*) (cosh 62 — cos p) (sin tp — i sinh x)) 
The non-symmetric parametrization given by Eqs. (|5.1| - |5.12 ) looks rather complicated. The 



calculation of the Jacobian is most lengthy but this has to be done only once. At the same 
time, the Jacobian does not contain Grassmann variables and the free energy F [Q] , Eq. 
( |5.13| ), is simple enough. Moreover, the supermatrix Q can be written as in the preceding 
Section in the form of Eq. ( |4.19|) (although the supermatrix Q is now different from that 
for the unitary ensemble). This allows to integrate first over the matrix u and obtain Eq. 

Further simplifications come from the fact that as previously one obtains an uncertainty 
of the type Ox 00 because the integrand in Eq. (|4.20|) does not contain the variables k, k* 
whereas the Jacobians Jg, Eq. (|4.17|) , and J M , Eq. (|5.10|) , are singular at 0,61,62, fi — > 0. 



We have seen in the preceding Section that the uncertainties can be rather easily avoided 
and, as a result, one obtains a more simple integral. The "regularization" procedure, Eqs. 
( |4.21|44.26| ), led to the integral, Eq. ( |4.27|) , that contained the variables <p and x only. 



Similar transformations can be performed for the orthogonal ensemble. Proceeding as 
for the unitary ensemble let us introduce the function F mn [Q] 

F mn = F [Q\ - mSTr (ATAT) - nSTr (t 3 Yt 3 y) (5. 14) 



The second term in Eq. ( |5.14| ) can also be written in the form of Eq. ( 4.24| ). Using Eqs. 
( |5.2| - |5.7| ) we can write the third term as 



— nSTr (r 3 rr 3 r) = An (coshfl 2 - cos/x) (5.15) 



In analogy with the transformation of the integrand in Eqs. ( 4.21 , [4.22 ), we can represent 
exp (— F [Q]) as follows 

e -F = e -F mn + e -F m0 ^ _ e -F(^ + e ~F Qn ^ _ e -FMj + -F ^ _ ^Nj ^ _ & - F ^ 

(5.16) 

where 

pH = F - Fn F (n) — F - F n 

- 1 - 1 mn - 1 On; - 1 - 1 van - 1 mO 

The parameters m and n in Eqs. ( |5.15| , |5.16|) are arbitrary. Therefore, substituting Eq. 
( p,16| ) into Eq. ( [4.3| ) we can take the limit m, n — > 00. The contribution coming from the 



last term in Eq. ( |5.16| ) vanishes because all singularities are compensated for any m and n 



but the integrand does not contain the anticommuting variables k,k*. The limit m — > 00 
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allows to expand the supermatrix T, Eq. (|4.23|) , near 1 (and the supermatrix Q near A). 



As has been explained in the preceding Section, in the limit m — > oo one can replace T — > 1 
everywhere in the integrand omitting simultaneously J^dnn* in the elementary volume [dQ]. 
The same is correct now and one should remove Jgdndn* from [dQ] , Eq. ( |5.9|) (changing the 
sign). 

The other singularity at 9 2 , [i — ► in the first and third terms in Eq. (|5.16|) can be 
avoided in a similar way. In the limit n — > oo the supermatrix Y, Eqs. ( |5.2| - |5.7| ), is also 
close to 1. To make an expansion in small deviations Y from 1 one can use the following 
parametrization 




Y= (\-iX)(l+iX)-\ X=\ f (5.17) 



The blocks A and L satisfy the constraints A = —A, L = L, {A, r 3 } = 0, {L, r 3 } = 0, where 
{...} is anticommutator. These blocks can be written in an explicit form as 





where the 2x2 matrices / and Zcontain conventional complex numbers / and I, whereas £ 
and ( consist of anticommuting variables £ and £. The explicit form of these matrices is 



/= 7 ,!= ,n ,C= (5.19) 






In Eq. (PD , / is an arbitrary complex number, while for / one should integrate over the 
domain Imf > 0. 

Substituting Eqs. ( FT^oTT^ ) into Eq. (g!3D one should expand the term STr ^3^X3!^ 



up to quadratic terms in X and replace Y by 1 everywhere else in the integrand. Calculating 
the Jacobian we can see that the factor J^dRi B dR 1F should be replaced by 1. Of course, 
this concerns only the first and the third terms in Eq. ( |5.16| ) because the second term does 
not lead to any singularity in the integrand at 82 = fi = 0. In fact, the contribution from 
the third term in Eq. ( |5.16 ) is zero because it is not singular at 9 = 61 = and does not 



contain the variables k, k*. At the same time we understand what to do with the singularity 
at 62 = \i = 0. 

The result of this discussion can be formulated finally as follows. We should replace Eq. 
(©by 

P(e,y)=pV(e,y) + pW(e,y), (5.20) 

(6, y) = ~ m lim M / A [Q] exp (-F mn [Q]) dQ (5.21) 

P {2) = -^Jiim^J A[Q] (exp (-F m0 ) - exp (-F mn )) dQ (5.22) 
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The integrand in Eq. ( 5.21]) has both singularities. Therefore, one has to replace every- 
where in the integrand T and Y by 1 simultaneously replacing J r gJ '^dudn* dRisdRip in the 
elementary volume [dQ], Eq. ( |5.9| ) by —1. As concerns Eq. ( |5.22| ) the integrand has only 
the singularity at 9 = 9\ = and one should replace by 1 the supermatrix T only. In the 
elementary volume Jedndn* should be replaced by —1. 

The subsequent manipulations are rather straightforward. Integrating over the super- 
matrix u one obtains for P^ (e,y) and (e, y) analogs of Eq. ( |4.2U| ). Then, the function 
P^ (e, y) is expressed in terms of the integral over the variables t = sin ip and z = sinh \ 



PW (e,y) 



AAdx 1 



-a 2 (t 2 +z 2 ^-x(t-iz) 



4t 2 dtdz 
(t 2 + z 2 f 



(5.23) 



In the integral in Eq. ( |5.22| ) one has to integrate first over the variables p, p*, a, and a* and 
then, the function P^ (e, y) reduces to 



P (2) (e,y) 



4A dx 2 



^-a 2 (t 2 +z 2 ) -x(Uj-i\z) 



(t-izfxH* 
(t 2 + x 2 ) 2 



dtdzdudX 



(5.24) 



where u = cosp, and A = cosh 62. The integration in Eq. ( |5.23| , p.24[ ) is performed over 
t and z in the intervals — 1 < t < 1, — 00 < z < 00 and over uj and A in the intervals 

-1 < LO < 1, 1 < A < OO. 

The integration over u and A in Eq. ( 5.24 ) can be carried out immediately. However, 
to provide the convergence of the integral over A one should shift the contour of integration 
over z into the complex plane z — > z + i5sgn (x), where 5 is an infinitesimal positive number 
and 

/ \ / 1, x > 
son (x) = < 1 „ 
y w [-1, x <0 

Intergrating over u and A and adding Eqs. ( 5.23| , 5.24 ) we obtain for P(e,y), Eq. ( |5.2C| ) 

v d 2 l{x) 



P(e,y) 



4A dx 2 



(5.25) 



I (x) = f 1 r e- a2 ( t2+z -) le-'W (t + ^_) 2 - e< t+i ^ (t - z^) 2 l (5.26) 

K J J-1J-00 L v ' v 1 lz _ (f + z lf v ; 

where ^_ = z + i5sgn (x). 

It is clear from the form of the function / (x) that it is convenient to differentiate first 
over x and then calculate the integral. However, one should be careful performing this, at 
first glance trivial, manipulation. The problem is that z_ contains x, which can result in an 
additional contribution. 

To avoid lengthy calculations let us consider first the case when x is finite nonzero 
number. Then, the derivatives dz-/dx and d 2 Z-/dx 2 vanish and one has to differentiate the 
exponentials only. Shifting the contour of integration z — > z+ffi, which can be done without 
crossing singularities in the complex plane and changing the new variable z as z — > z/a one 
obtains 
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V 



Pc (e, y) = — exp 



x 



Aa? 



J xt sinh xt exp a 2 t 2 ) 



00 exp (— z 2 ) 



2 1 ^1 

Z + 4a 2 



(5.27) 



(the variables 2 and ?/ are related through Eq. ( |4.2j )). 

Eq. (|5.27| ) holds for any finite x but is it the final result? It would be the final result 
it the density function were continuous at x — 0. As concerns the unitary ensemble, we 
already know that the function P (e, y)is continuous (see Eq. (|4.28 )) but does the continuity 
follow from a physical principle? In fact it does not and the function P (e, y) for the unitary 
ensemble contains a 5-function at x — 0. 

To extract the 5-function let us expand the exponentials in the integrand in Eq. G5.26|) . 
In the first two orders one obtains 



P(e,y) 



v_£_ 
2Adx 2 



t 2 e 



-a 2 (t 2 +z 2 _) 



X 



{t 2 + zlf iz-t 2 + z 2 __ 



dtdz 



(5.28) 



The first term in the integrand in Eq. ( |5.28| ) has no singularities and one can shift the 
contour of the integration over z such the variables z_ are replaced by z. Then, this part 
of the integrand does not contain x and the differentiation gives zero. The contribution 
involved comes from the second term in the integrand. Writing zZ 1 as 

1 z — iSsgn (x) 



z_ z 2 + 5 2 

one can represent the function P (e, y) for x — > as 

V d 2 rl rco 



x-*0 



a h m 

2A dx 2 s^o 



a 2 (t' 2 +z 2 ) t \X\ 
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t 2 + z 2 z 2 + 5 2 



dtdz 



The integration over z in the limit 5 — > is elementary and one obtains for the anomalous 
contribution P r (e, y) the following expression 



P r (e,y) 



Itxv 



S (x) J exp (— a 2 t 2 ^j dt 



(5.29) 



Making some simple transformations in Eq. ( [5.27] ) the final result for the density of complex 
eigenvalues P (e, y) can be written as 



P(e,y)=P r (e,y) + P c (e,y) 
where P r (e, y) is given by Eq. ( |5.29[) and P c (e, y) equals 

2ttu 1 



P c (e,y) = r-^j £ t sinh (\x\t) exp (-a¥ 



dt 



(5.30) 



(5.31) 



where $ (v) = ^= exp (— u 2 ) du. It is not difficult to check that the function P (e, y), Eqs. 
( |5.29| - |5.3l| ), satisfies the normalization condition, Eq. ( |4.29| ), and the singular part P r (e,y) 
gives an essential contribution that becomes small only in the limit a — » 00. The function 
P c (e, y) is represented in Fig.0. 
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The existence of the anomalous part P r (e,y), Eq. ( |4.32|) , means that a finite fraction of 
all eigenvalues remains real for any imaginary vector potential h in the models of disorder, 
Eqs. (JO, p.3| ) or degree of asymmetry a for the real random matrix models. At the same 
time, the function P c (e, y) decays when y — > 0, which corresponds to a vanishing probability 
of eigenstates with small but nonzero imaginary parts. 

In contrast to the unitary ensemble, the function P (e, y) for a <C 1 can hardly be 
obtained from a perturbation theory. Most of the eigenvalues are in this case real. In the 
opposite limit a ^> 1 one should distinguish between several regions. In the limit \x\ <C a 
the asymptotics is determined by the expression 

P c e, y) ~ - - v 11 5.32 
v ,yj 2a 2 A 2a K ' 

showing a linear decay of the density as \x\ — > oo. 

In the region |x| 3> 2a the density of complex eigenvalues is constant for \x\ < 2a 2 and 
falls off outside this interval. Its value in this region is the same as in the unitary case, 
Eq. (|4.31|) . This corresponds to the elliptic law, Eq. (|4.32|) . For an ensemble of strongly 
asymmetric real random matrices with a Gaussian distribution this law has been proven in 
Ref.0H. The authors of this publication have also found numerically that the portion of 
real eigenvalues for their ensemble decays as iV~ 1//2 , where N is the size of the matrices. 
Apparently, this behavior corresponds to the ^-functional part P r (e,y), Eq. ( |5.29| ), in the 



eigenvalue density for the case of weak asymmetry (orthogonal analog of Eqs. ( |4.6| , |4.6|)). 



VI. DISCUSSION 

The results presented in the previous Sections demonstrate that the disorder models with 
a direction are interesting and can be efficiently studied using the supersymmetry technique. 
The a- model derived, Eq. ( |3.29| ), can be used in any dimension. It is relevant to emphasize 
that, as usualSi, the dimensionality is determined by the geometry of the sample. So, the 
one-dimensional version of the a-model corresponds to a thick wire with a directed hopping. 
In the language of vortices in a superconductorlll the ID model can describe the vortices 
in a slab with line defects and the magnetic field parallel to the surface. Such a model is 
somewhat more realistic than a purely ID model of Ref.0. The 2D er-model is supposed 
to describe the vortices in a bulk superconductor with line defects. In addition, one can 
imagine a situation when the sample is long but has a small cross-section. If the line defects 
are aligned in the longitudinal direction one comes to the 0D a-model considered in the 
present paper. 

Of course, the directed non-Hermitian Hamiltonians can arise not only from the vortex 
model but also correspond to non-equilibrium processes. A very interesting possibility is the 
directed hopping model, Eq. (|2.3j ) that can be considered as a quantum counterpart of the 
directed percolation modeEiJ. Applications to other physical systems that can be reduced to 
models of a disorder with a direction also deserve an attention. The problem of turbulence 
is one of most famous. The main features of the turbulence are believed to be described by 
the Burgers equationH'H'0. Reduction of the Burgers equation to a linear equation allows 
to use well developed methods of disorder physics. A similarity of the linear equation to 
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equations used in study of problems of directed polymers have already inspired application 
of the replica method to study the problem of turbulencel!. Use of the supersymmetry for 
the problems of the turbulence might be one more interesting direction of research. 

Leaving these interesting problems for future study let us summarize the results obtained 
in the present work. The cr-model, Eq. ( |3.29| ), differs from the a-models used in the localiza- 
tion and mesoscopic problemsi'i by the term with the matrix Ai. Although the Hamiltonians 
with the direction, Eqs. ( |2.1| , |2.3| ) can be obtained from conventional Hermitian Hamilto- 
nians in a magnetic field by the formal replacement A — nh, the same replacement in the 
conventional cr-models would not lead to Eq. ( |3.29j ). This reflects an essential symmetry 
difference between systems in a magnetic field where the time reversal invariance is broken 
and the models with direction that are time reversal invariant. 

In contrast to average density of states for Hermitian disorder problems which is always 
smooth, the joint probability density of complex eigenenergies considered in the previous 
Sections is a non-trivial quantity. The cr-model was derived to describe this quantity and 
it is expected to be sensitive to localization-delocalization transitions in one- and higher 
dimensional systems^. 

The the form of OD version of the cr-model obtained above demonstrates the equivalence 
between the directed disorder models in a limited volume and ensembles of random weakly 
non-Hermitian or weakly asymmetric real matrices that have been mapped onto the OD cr- 
model previously]^. Complex random non-Hermitian matrices appear in study of dissipative 
quantum maps@'0 whereas random real asymmetric matrices have applications in neural 
network dynamicJiH. So, the cr-model can describe completely different phenomena in an 
unified manner. 

The supermatrix cr-model can serve as an useful calculational tool for all these non- 
Hermiatian problems. Although the new term with the matrix Ai in the cr-model, Eq. ( ^.29[ ), 
makes the use of previous parametrizationsi difficult, the new parametrization suggested in 
the present paper allows to circumvent the difficulties and obtain in a straightforward manner 
explicit results for the OD case. Weakly non-Hermitian random matrices can also be studied 
using more traditional methods of orthogonal polynomials^. However, study of weakly 
non-symmetric real matrices with this method seems be more difficult and the density of 
complex eigenvalues, Eqs. (|5.29| - |5.3ll) , has been calculated for the first time. Besides, the 
cr-model approach is not dependent on details of the model considered and can be applied 
not only to Gaussian models. It can also be used to study the directed models in one and 
higher dimensions where one can expect localization-delocalization transitions. 

Eqs. ( |5.29[ - p3ID demonstrate that at any finite disorder and "imaginary vector-potential" 
a finite portion of eigenvalues remain real whereas this does not occur if the time reversal 
invariance is broken, Eq. Q4.28 ). This phenomenon has manifested itself in numerical study 
of different models. In Reiki! ^1 ensembles of random strongly asymmetric matrices (sym- 
metric and antisymmetric parts had the same order of magnitude) were considered. It was 
found that the fraction of real eigenvalues decayed as N~~ l l 2 for large matrix sizes N. Ap- 
parently, this corresponds to the finite fraction of the real eigenvalues P r (e,y), Eq. ( |5.29| ). 
because in the ensemble of weakly non-symmetric matrices involved the magnitude of the 
antisymmetric part of the random matrices is N 1 / 2 times smaller than that of the symmetric 
one. 

A finite fraction of real eigenenergies was found in a numerical study of the 2D model, 
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Eq. ( |2.3| ), (without magnetic interactions) near the center of the banofi Although the 2D 
case was not considered in the present paper and nothing can be said about a possibility of 
a mixture of eigenstates with real and complex eigenvalues one can argue that, may be, the 
parameters of the model of Ref.El corresponded to the OD case. This might easily happen 
because the localization length in weakly disordered 2D systems is exponentially large and 
can exceed the sample size, which would correspond to the OD regime. If this is really so, 
the results of the present study are in an agreement with the numerical investigation. 

The phenomenon that some finite portion of eigenvalues lies on a certain line in the 
complex plane occurs also in other models with a randomness. Recently, it was found that 
a finite fraction of all roots of random self-inversive polynomials lies on the unit circled. 
At the same time, if the polynomials are not self-inversive the density of complex roots in 
smooth everywhere in the complex plane. 

It is clear from the preceding discussion that the directed disorder models deserve further 
investigation. 

VII. APPENDIX 

A. Non-Hermitian parametrization (unitary ensemble) 

Let us calculate for the unitary ensemble the Jacobian for the parametrization given by 
Eqs.. (f3~7 - 4.1lT ) (it was suggested to call it "non- Hermit ian parametrization". As usuaffl, 



it is convenient to consider the length Str (dQ) 2 . With Eq. ( f4.7| ), it can be written as 

STr (dQf = STr ((dQ ) 2 + [5T, Q f + ^T5Q ) (7.1) 

where 5T = TdT, 5Q = QodQ and [., .] is the commutator. 
It is easy to see from Eq. (f4.8|) that 



and hence 



{5Q , A 1 } = (7.3) 



where {., .} is the anticommutator. 

Then, using the relation [ST, Ai] = and Eq. ( [7.3| ) we obtain 

STr (5T5Qo) = STr (A^TSQaAj = -STr (5T5Q ) = (7.4) 

which shows that Jacobians is the product of Jacobians corresponding to dQo and 5T. As 
concerns dQo, we have 



STr(dQ ) 2 = A((d V ) 2 + (d X ) 2 ) (7.5) 



Writing Eq. (g^) as 
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T = uT v 



(7.6) 



one obtains 



ST = vT SuT v + v5T v + Sv 



(7.7) 



where, with Eq. ( f4.10| , [4.11 



5Tn = - 



i_ ( d9 
2 { d§ 



(7.8) 



Su = 8u\\ + Su±, 



Su\\ — 2r 3 (rjdrj* — drjr]*) , Suj 



dr] 
-dfj 



and similar equations can be written for Sv. 

Substituting Eq. Q4.10] , |4.11| , |7.8|) into Eq. ( |7.7|) one can represent the supermatrix ST 

as 



ST = ST" + ST A 



(7.9) 



I, g - igi / d?/\ / d K 

5T "=2co8^-( _ d _ 1+2 1 _ d _ Q 



(7.10) 



+2r3 (77CZ77* — drjrj* + «d/«* — dnn*) + 4r 3 cos (K*dr/ — drfn) 



(7.11) 



1 / d0(l -4kk* 

2 I 



•j/i /1 , /i *\ +4 sin — - — (k di] + di] Kj) 
idQ\ (1 + 4kk*) / 2 



In Eqs. ([779- 7. 11 ) 5T" commutes with A, ST 1 - anticommutes with A. The second line in Eq. 
( [7. 1C| ) does not contribute to [ST, Qq] in Eq. (J77J). In Eqs. ( [7.1 U| , |7.11| ), one can change the 
variables 



d# (1 - Akk*) -> d0, d0i (1 + 4kk*) -> cZ6»i 



(7.12) 



and make the shifts 



i(de o \ i(de 

2 [ zd^ J ~* 2 I 



+ 4 (K*drj + dr]*K) sin (7.13) 



2 sin dr/ — > 2 sin dr/ + (d9 — id9%) k 



(7.14) 
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d ^ dK _ cos ^l d , (7 . 15) 



The transformations, Eqs. ( |7.12 - 7TT5| ) do not change the Jacobian and 5T" and 5T 1 - take 
a more simple form 



where zr3<ic is the second line of Eq. (|7.10|) and 1 is the unit 8x8 matrix. 
Further computation is already simple. Changing once more 

sin -dr] — > drj (7-18) 



one obtains a contribution to the Jacobian proportional to Jg, Eq. ( |4.17|) . Writing in the 
new variables the second term in Eq. ( |7.1| ) we have 

STr [5T, Q } 2 = 4 ((dO) 2 cos 2 <p + (c^) 2 cosh 2 X ) (7.19) 

+128 (cos 2 ^^drjdv* + sin 2 — — ^— — dhidK* \ 
Eqs. ( |7.5| , |7.19|) lead to the elementary volume [dQ], Eq. ( |4.15| ). 



B. Non-symmetric parametrization (orthogonal ensemble) 

To calculate the Jacobian of the parametrization, Eqs. ( |5.1[ - p.7| ), for the orthogonal 
ensemble we can use the results obtained for the unitary ensemble because Eq. ( |5.1|) contains 
the same supermatrices Qo and T as previously. However, the presence of the supermatrix Y 
makes the computation quite lengthy. The length STr {dQ) 2 is written as in the preceding 
Subsection 

STr (dQ) 2 = STr ((dQ ) 2 + [8Z, Q } 2 + A5Z5Q ) (7.20) 
where 5Z = ZdZ can be written as 

5Z = SR (Y 5TYo + 5Y + dRR + RdSSR) RS (7.21) 



The last term in Eq. (|7.20| ) is equal to zero (see Eq. (|7.4|) ). As concerns the supermatrix 



<5T, it can be written after the replacements, Eqs. ( |7.12| - |7TT5| ), in the form of Eqs. ( |7.16| . 
~17|) . So, one has to calculate the other differentials entering Eq. ( |7.21| ). Using Eq. 



one can rewrite dSS to the form 
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dSS = (dSS) + (dSS) ± , (7.22) 

(dSS) = 2r 3 l (daa* - ada*) , (dSS) ± = 2iA 1 da 
Taking the supermatrix R from Eq. ( |5.6|) one can derive 

RdSSR = R(dSS) ± R+(dSS) (7.23) 

R (dSS) ± R = 2%k x da + 4iAi (dap* + pda*) (7.24) 

and 

dRR = 1 (2 ( _° _ q P J + 2r 3 (dpp* - pdp*)^ (7.25) 

Now we have to calculate 5Y Q . Using Eqs. ( |5.3| - |5.5| ) one can represent this differential in the 
form 

5Y = 8Y 1 + 5Y 2 + Y 1 Y 2 SY 3 Y 2 Y 1 (7.26) 
Calculating the matrices SYi, SY 2 , and 8Y3 we rewrite SYq as follows 

5Y = l(U d ^ W AR , Q )--dph °Y) 

\2 \ rf/?ir 3 cosh^ 2 J 2 J J 



where 



/ — i \ ( cos u — sin u 

t 2 = . , wr 3 w = 

u I \ — sin /x — cos fi 



Making the replacement 

and the same for dry and dr?* one can derive 

Y Q 5TY = 21 (cosh |( J k , f ) + , smh I ( ^ ) ) ,7,8, 

o • a ( e * ( drf \ . . , 2 ( dn'n \ i A 

+2tAi cos — ,_. ' — zsinn — ,_, „ dv 

\ 2 \ drf j 2 \ nek' j 2 / 

where dr/ = weir?, ci??' = dfjw and the same for dn' and d«. The contribution from ir 2 dc, Eq. 
( [7.161 ), is not written because it can be removed by a proper shift of d/3 and dp. 
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Substituting Eqs. (|7.22[ - |7.25|) into Eqs. ( |7.21| , |7.20|) we see that second terms of Eq. 
( 7.25 , |7.25|) do not contribute. After making the replacement in ST, Eq. ( 7.18 ), and shifting 



da = dui — (cosh — cos —dri + i sinh — sin — dn* 
\ 2 2 1 2 2 



(7.29) 



dp = dpi — ( cosh cos ^g?k — i sinh sin ^drf 



it is convenient to introduce the matrix differentials 
, / da i da 2 \ , / 

da=[ -da; -dot 



dpi dp 2 
-dp* 2 —dpi 



where 



dao 



i ^2 . ft j * • • 1^2 ft, 

— cosh — sin —dri — i sinh — cos —an 
2 2 1 2 2 

r * ■, &2 ■ ft , ■ . -, &2 ft , * 

da = cosh — sin —dri — % sinh — cos — an 
2 2 2' 22 

^2 ft 62 ft 

dp 2 = — cosh — sin — an + 1 sinh — cos — dr\ 

dp% = cosh — sin —dn + i sinh — cos — dri* 
H2 2 2 2 2 1 



The Jacobian J M of the transformation, Eqs. (|3.17|) , equals 

4 



Then Eq. ( |7.21| ) can be written as 



(cosh 6 2 — cos nY 



(7.30) 



(7.31) 



(7.32) 



5Z = SR5URS, 



(7.33) 



SU = 5Y + zAi ( 2do- - -d(9 ) + 2A;dpl (7.34) 



2 



where the matrices da and dp entering da and dp, Eq. (|5.7| ) have the structure Eq. (|7.30|) 
and 

'i-i) 



fc 



. One can calculate 5Z, Eq. ( |7.33|) , calculating first R5UR and then SZ. The corresponding 
manipulations are still quite lengthy. One should again make different replacements that 
do not change the Jacobian. Alternatively, one might write the final result using general 
symmetry properties of 5Z. Finally one obtains 

5Z = 5Y^ + iAx (2da - -d§) + 2kdpl (7.35) 
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The supermatrix SYq entering Eq. (|7.35|) equals 



+ dp 



sinh (9 2 c?/3i 





r 2 

l2 



de 2 



r 2 0\ 
0/ 

f 

7i 



(7.36) 



Using Eq. ( |7.35| ) we can calculate STr [5Z, Q ] . The contribution of the anticommuting 
part SZ± decouples from that of the commuting one 5Z\\ and one obtains 



STr [5Z±, Qo] = Q4[da\dal (1 + cos (p + £%)) + da 2 da 2 (1 + cos (p — £%)) 

+dpidp{ (1 — cos (<p — z%)) + dp 2 dp* 2 (1 - cos (<p + 
The Jacobian J ra corresponding to the length, Eq. (|7.37|) equals 

1 1 



(7.37) 



J, 



024 / . o . , o \ z 

( sin <y9 + sum x J 
The commuting part 5Zn contributes to the elementary length as 



(7.38) 



STr 



SZ\\ , Qo = 4[((tZ/if + (d,3) sin' //) sin' 



(7.39) 



+ (d#) 2 cos 2 9? + {ddxf cosh 2 x + {d6 2 f + (d/^) 2 sinh 2 2 ] 

Combining the contribution coming to the Jacobian from Eqs. (|7.5| , [7. 39| ) with those 
written in Eqs. ( |7.32| , |7.38| ) and recalling that the replacement, Eq. ( |7.18| ), results in an 
additional multiplier proportional to Jg one obtains finally the elementary volume [dQ], Eqs. 



34 



REFERENCES 



1 P.W. Anderson, Phys. Rev. 109, 1492 (1958) 

2 B.L. Altshuler, Sov. Phys. JETP Lett. 41, 648 

3 PA. Lee and A.D. Stone, Phys. Rev. Lett. 55, 1622 (1985) 
4 R.G. Laughlin, Phys. Rev. B 23, 5632 (1981) 

5 L.P. Gorkov, A.I. Larkin, and D.E. Khmelnitskii, Sov. Phys. JETP Letters, 30, 228 (1979) 
6 E. Abrahams, P.W. Anderson, and T. V. Ramakrishnan, Phil. Mag. 42, 827; E. Abrahams 

and T.V. Ramakrishnan, J. Non-Cryst. Sol. 35, 15 
7 K.B. Efetov, Adv. in Phys. 32, 53 (1983) 

8 K.B. Efetov, Supersymmetry in Disorder and Chaos, Cambridge University Press, New 
York (1997) 

9 J. J.M. Verbaarschot, H.A. Weidenmiiller, and M.R. Zirnbauer, Phys. Rep. 129, 367 (1985) 

10 O. Bohigas and M.-J.Gianonni, in Mathematical and Computational Methods in Nuclear 
Physics, ed. J.S. Dehesa, J.M.G. Gomez, A. Polls, Lecture Notes in Physics, Vol. 209, 
Springer- Verlag, Heidelberg (1984) 

11 O. Bohigas, in Chaos and Quantum Physics, ed. M.-J. Gianonni, A. Voros and J. Zinn- 
Justin, Elsevier, Amsterdam (1991) 

12 G. Casati, ed., Chaotic Behavior in Quantum Systems, Plenum Press, New York (1985) 
13 M.L. Mehta, Random Matrices, Academic Press, San Diego (1991) 

14 F. Haake, Quantum Signatures of Chaos, Springer, Berlin (1992) 
15 B.A. Muzykantskii and D.E. Khmelnitskii, JETP Lett. 62, 76 (1995) 

16 A. V. Andreev, O. Agam, B.D. Simons, and B.L. Altshuler, Phys. Rev. Lett. 76, 3947 
(1996) 

17 M.R. Zirnbauer, Nucl. Phys. A560, 95 (1993) 

18 N. Hatano and D. Nelson, Phys. Rev. Lett. 77, 570 (1996) 

19 L. Civale, A.D. Marwick, T.K. Worthington, M.A. Kirk, J.R. Thompson, L. Krusin- 
Elbaum, Y. Sum, J.R. Clemm, and F. Holtzberg, Phys. Rev. Lett. 67, 648 (1991); M. 
Konczykowski, F. Rullier-Alebenque, E.R. Yacoby, A. Shaulov, Y. Yeshurin, and P. Lejay, 
Phys. Rev. B44, 7167 (1991); R.C. Budhani, M. Swenaga, and S.H. Liou, Phys. Rev. Lett. 
69, 3816 (1992) 

20 D.R. Nelson and V. Vinokur, Phys. Rev. B48, 13060 (1993), and references therein. 

21 S.P Obukhov, Physica, A 101, 145 (1980) 

22 J.M. Burgers, The Non-Linear Diffusion Equation, Reidel, Dordrecht (1974) 

23 A. Polyakov, Nucl. Phys. B396, 367 (1993); Phys. Rev. E 52, 6183 (1995) 

24 M. Kardar, G. Parisi, and Y.C. Zhang, Phys. Rev. Lett. 56, 889 (1986) 

25 J.P Bouchaud, M. Mezard, and G. Parisi, Phys. Rev. E 52, 3656 (1995) 

26 H.C. Fogedby, A.B. Eriksson, and L.V. Mikheev, Phys. Rev. Lett., 75, 1883 (1995) 

27 V. Gurarie and A.A. Migdal, Phys. Rev. E 54, 4908 (1996); G. Falkovich, I. Kolokolov, 
V. Lebedev, and A.A. Migdal, Phys. Rev. E 54, 4896 (1996); E. Balkovsky, G. Falkovich, 
I. Kolokolov, and V. Lebedev, Phys. Rev. Lett. 78, 1452 

28 J. Ginibre, J. Math. Phys. 6, 440 

29 V.L. Girko, Theor. Prob. Appl. 29, 694 (1985); Theor. Prob. Appl. 30, 677 (1986) 
30 R. Grobe, F. Haake, and H.-J. Sommers, Phys. Rev. Lett. 61, 1899 (1988); R. Grobe and 
F. Haake, Phys. Rev. Lett., 62, 2893 (1989) 



35 



31 H.-J. Sommers, A. Crisanti, H. Sompolinsky, and Y. Stein, Phys. Rev. Lett. 60, 1895 
(1988) 

32 N. Lehmann and H.-J. Sommers, Phys. Rev. Lett. 67, 941 (1991) 

33 H. Sompolinsky, A. Crisanti, and H.-J. Sommers, Phys. Rev. Lett. 61, 259 (1988) 

34 B. Doyon, B. Cessac, M. Quoy, and M.Samuelidis, Int. J. Bifurc. Chaos 3, 279 

35 R.A. Janik, M.A. Nowak, G. Papp, and I. Zahed, preprints |cond-mat /96 1 2249| , |hep- 



ph/9606329|; R.A. Janik, M.A. Nowak, G. Papp, J. Wambach, and I. Zahed, preprint 



[hep-ph/9609491 
36 J. Feinberg and A. Zee, preprint hep-th/9703087 

37 Y.V. Fyodorov, B.A. Khoruzhenko, and H.-J. Sommers, Phys. Lett. A 226, 46 (1997) 

38 K.B. Efetov, preprint, |cond-mat/ 970209 1| 

39 A. Altland, S. Iida, and K.B. Efetov, J. Phys. A 26, 3545 (1993) 

40 A. Chekhlov and V. Yakhot, Phys. Rev. E 51, R2739 (1995) 

41 Y.V. Fyodorov, B.A. Khoruzhenko, and H.-J. Sommers, |cond-mat/9703152| 

42 E. Bogomolny, O. Bohigas, and P. Leboeuf, J. Stat. Phys. 85, 639 (1996) 



36 



FIGURES 



FIG. 1. The density of complex eigenenergies P (e, y) for the unitary ensemble as a function of 
the imaginary part (x = 27ry/A) for a = 1, 2, 3 



FIG. 2. The density of complex eigenenergies P c (e, y)for the orthogonal ensemble as a function 
of the imaginary part (x = 27ry/A) for a = 3, 5, 7 
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